Sign changes and resonance of intrinsic spin Hall effect in two-dimensional hole gas 
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The intrinsic spin Hall conductance shows rich sign changes by applying a perpendicular mag- 
netic field in a two-dimensional hole gas. Especially, a notable sign changes can be achieved by 
adjusting the characteristic length of the Rashba coupling and hole density at moderate magnetic 
fields. This sign issue may be easily realized in experiments. The oscillations of the intrinsic spin 
Hall conductance as a function of 1/B is nothing else but Shubnikov-de Haas oscillations, and the 
additional beatings can be quantitatively related to the value of the spin-orbit coupling parameter. 
The Zeeman splitting is too small to introduce effective degeneracy between different Landau levels 
in a two-dimensional hole gas, and the resonant intrinsic spin Hall conductance appears in high hole 
density and strong magnetic field due to the transition between mostly spin- — | holes and spin-| 
holes is confirmed. Two likely ways to establish intrinsic spin Hall effect in experiments and a 
possible application are suggested. 

PACS numbers: 72.25.-b,72.25.Dc 



In the developing field of spintronics 1 -, which is be- 
lieved to be a promising candidate for future informa- 
tion technology^, the generation and efficient control of 
spin current is a key issue. So the intrinsic spin Hall ef- 
fect (ISHE) predicted by Murakami et al£ and Sinova et 
al£ has generated intensive theoretical studies and this 
ISHE is associated with the dc-ficld-induced transitions 
between the spin-orbit-coupled (SOC) bands, and con- 
tributes from all occupied electron states below the Fermi 
energy. The ISHE is distinguished from the extrinsic spin 
Hall effect (ESHE), which is due to impurity scattering 5 . 
On the experimental side, the spin Hall effect (SHE) re- 
ported in a two-dimensional hole gas (2DHG) is likely of 
the intrinsic origin^; nevertheless, the intrinsic or the 
extrinsic origin of the SHE observed in a 3D electron 
film& is still under dabat o 9 ' 10 . 

The intrinsic spin Hall conductance (ISHC) of a 2DHG 
in a perpendicular magnetic field has been studied within 
fc-cubic Rashba model by M. Zarea et oJJ-L and a general- 
ized Luttinger model by ua— most recently. The work by 
M. Zarea et al. only focused on the low field regime, and 
a remark 1 ^ has been given by J. Schliemann on this as- 
pect: in any case, influence of a magnetic field coupling to 
the orbital degrees of freedom should only be appreciable 
if the field is strong enough to produce typical cyclotron 
radii being of order of the system size or smaller; there- 
fore, arbitrarily small fields cannot be expected to have 
an effect in real experiments. Within a more general band 
structure Hamiltonian, we have predicted that the "effec- 
tive" energy crossing near Fermi energy between mostly 
spin- — i holes and spin-| holes at a typical magnetic field 
gave rise to a resonant ISHC™ with a certain hole den- 
sity, however, a strong magnetic field region is required. 
Except resonanc e) 12 ' 14 : 15 , the sign issue is another crucial 



aspect of ISHGiSiii, which can be modulating by exter- 
nal magnetic fielcUS. To be convenient for future exper- 
imental detection or possible applications, it is of course 
desirable to discuss the effect of magnetic field in detail, 
especially in moderate magnetic field region, which can 
be easily achieved in experiments or technologies. 

In this paper, we will concentrate on the sensitivity of 
ISHC to magnetic field, the length scale of the Rashba 
coupling, and hole density, which shows rich sign changes. 
We predict that a more notable sign changes can be 
achieved by adjusting the length scale of the Rashba 
coupling at certain moderate magnetic fields for a fixed 
hole density. By discussing the effect of Zeeman split- 
ting on energy levels, we clarify that the interplay be- 
tween mostly spin — | holes and spin-| holes is the ori- 
gin of resonant ISHC in a 2DHG under a perpendicu- 
lar magnetic field, which is different from the case 1 ^ in 
two-dimensional electron gas (2DEG). In addition, the 
oscillations of the intrinsic spin Hall conductance as a 
function of 1/B is nothing else but Shubnikov-de Haas 
(SdH) oscillations, and the additional beatings can be 
quantitatively related to the value of the SOC parame- 
ter. Two likely ways to establish ISHE in experiments 
and a possible application are suggested, which may be 
easily realized in experiments. 

When the quasi 2DHG is sufficiently narrow and the 
density and the temperature are not too high, only the 
lowest heavy-hole subband is occupied. The effective 
Hamiltonian for a single heavy hole subjected a spin- 
orbit interaction due to structural-inversion asymmetry 
(SIA) can be written aa 18 : 19 : 20 
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using the notations p± = p x ± ip y , a± = (a x ± ia y )/2, 
where p, a denote the hole momentum operator and Pauli 
matrices respectively. These Pauli matrices operate on 
the total angular momentum states with spin projection 
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±3/2 along the growth direction, and we simply call it 
holes with spin-±| below. In the above equation, m is 
the heavy hole mass, and a is Rashba SOC coefficient 
due to SIA. The validity of above modeli£ii2i2S given by 
the Hamiltonian (1) was restricted to sufficiently small 
wave numbers and densities. 

We follow the method by M. Zareaii>ei al. and uaH 
bellow. We impose a magnetic field B = Bz by choosing 
the Landau gauge A = —yBx, and then p x = fak x + —y, 
where — e is the electric charge. The destruction operator 
a = , 1 (p x + ip v ), and the creation operator = 

i 1 - (Px—iPv) are introduced to describe Landau levels, 

V2mnuj x y ' 

where w = — . These operators have the commutation 
[a, a)\ = 1. In terms of these operators, the Hamiltonian 
in the absence of Zeeman coupling is 

H = hu[(a)a + -) + i 7 (a 3 o- + - a+V)], (2) 

in which the dimensionless parameter 7 is defined by 
7 = 2^y^^. In units of fau, the eigensystem can 
be expressed as 



— I COS I 

sin 6 



ns<Pn-3 



E ns = (n - 1) + sW 7 2 n(n - l)(n - 2) + -, (3) 

with s = ±1 for n > 3, s = 1 for n < 3 and <f> n is 
the eigenstate of the nth Landau level in the absence 

of SOC. Besides, 9 ns = arctanl 



27^/n(n-l)(n-2) 



l s 1 /47 2 n(n-l)(n-2)+9+3 J ' 

where ns G (0, for n > 3, and 6 = f when n < 3. 

When the system is driven by a dc electric field in — y 
direction, the most interesting case is the z-direction- 
polarized spin current along x direction, and the cor- 
responding spin current operator can be defineo»2*A£ as 
is = ^fr{vx,& z } where v x is the velocity operator. We 
shall only focus on the first order in the perturbation 
term because < n, s\js\n, s >= 0, which is 

usu = y (n> s| is K; s>) (n/> s/| w k s) + H.c, (4) 



n' s' 



Ens En's' 



where (js)ns is the current carried by a hole in the state 
|n, s > of Hq, and e ns = fajE ns . In addition, n' = n ± 1 
since the matrix element vanishes for other values of n' . 
The average current density of the TV^-hole system is 
given by I s = S„ s 0's)ns/(ens) where / is the Fermi 
distribution function, and Nh = J2 ns f( e ns)- The ISHC 
is then given by Gs = Is/EL X , where the Landau de- 
generacy factor eB/(hc) is also included. We can finally 
obtain the complete ISHC, and the 2DHG system we con- 
sider is confined in the x — y plane of an area L x x L y . 

Experimental study^ has been performed for heavy- 
hole density nj = 2 x 10 16 m~ 2 , and the characteristic 
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FIG. 1: (a) The ISHC as a function of 1/B for n h = 2 x 
10 16 m -2 , and the characteristic length am/h 2 = 0.25nm of 
the Rashba coupling 19 ", (b) Same as Figl.(a) but am/h 2 = 
0.075nm. 



length am/fa 2 = 0.25nm of the Rashba coupling—. By 
using these parameters, the ISHC as a function of 1/B 
is shown in Fig.l (a). The ISHC oscillates as a result 
of the alternative occupation of mostly spin-| holes and 
■| holes, which is consistent with our pre- 



mostly spin- 

vious result. The oscillations of the ISHC as a func- 
tion of 1/B are nothing else but another manifestation 
of SdH oscillations. Following the same argument of 
usual SdH2£, the periodicity of density of states with 
1/B results from the equality of the area of two suc- 
cessive electron orbits in fc-space and the oscillation fre- 
quency is given by f= nh ^° =41.36, where 4>o=- is the 
flux quantum. In relative weak magnetic field region 
{1/B > 0.22/T in Fig.l (a)), there are two indepen- 
dent oscillation frequencies originating from the two spin- 



splitting subbands, 
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^(kf) 2 , where n+(n h ) 
are respectively the hole density in spin up (down) sub- 
bands and krj are respectively their Fermi momenta. It 
is just the composition of these two frequencies produces 
the additional beating pattern appearing in our numer- 
ical result, and the additional beatings can be quantita- 
tively related to the value of the SOC parameter 

kf = - v 7 ! - 16nn h L 2 ) 

+ ( J - V 7 ! ~ l^n h L 2 ) + Snrih} 1 / 2 , (5) 

where L = am/h 2 . Fast Fourier transformation has been 
carried out separately on the strong- and weak magnetic 
field data shown in Fig.l and frequencies obtained agree 
numerically with Eq. (5) and f± = ^^(k^) 2 rather well. 

The importance of sign changes of ISHC has been em- 
phasize by Yao et alm& by applying first-principles calcu- 
lations to study SHE in semiconductors and simple met- 
als in the absence of a magnetic field. For the definition of 
sign, positive spin Hall conductance means that spin-up 
component flows to the positive x direction. The sign of 
ESHE induced in the skew scattering, which dominates 
over the side jump process in the weak disorder limit, 
depends on the sign of the impurity potential, and does 
not change with the carrier density or the Rashba cou- 
pling strength. As shown in Fig. 1(a), the sign of ISHC 
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FIG. 2: The optimal am/h 2 for the most notable sign changes 
as a function of hole densities (solid line) , and the correspond- 
ing most favorable range of magnetic field is also plotted as a 
function of hole densities (shadow). 



changes periodically between v — 6 and v = 13. Fur- 
thermore, this sign changes is robust in this case, and 
a more notable sign changes is sensitive to am/h 2 for a 
fixed hole density, which can be achieved by modulat- 
ing am/h 2 at certain magnetic fields. Fig.l (b) accounts 
for this. A more notable sign changes is shown for the 
same hole density as in Fig. 1(a) at am/h 2 = 0.075nm, 
where the magnetic field is in the range 2.07T ~ 3.60T, 
especially when 2.50T < B < 2.85T. 

To be convenient for future experimental detection, we 
have made an scan on the magnetic field, and the charac- 
teristic length scale dependence on the most notable sign 
changes as a function of hole densities, and the optimal 
am/h 2 for the most notable sign changes as a function of 
a hole density has been shown in Fig. 2, beside, the cor- 
responding most favorable range of magnetic field is also 
plotted as a function of hole densities, which denoted by 
shadow in Fig. 2. Since the ISHC in +x direction is larger 
than the ISHC in — x direction for most case, the opti- 
mal am/h 2 and the corresponding most favorable range 
of B for the most notable sign changes means that when 
the ISHC in —a; direction is the largest for a fixed hole 
density. 

Moreover, The calculated ISHC shows sign changes as 
hole density varies for a fixed magnetic field and the char- 
acteristic length scale of Rashba coupling, which differs 
from the ESHE too. As shown in Fig. 3(a), the impor- 
tance of this aspect is that, the required magnetic field 
(B = 0.827T) is easy to achieve in experiments, and the 
sign changes of ISHC may be detected expediently. It's 
a likely way to establish the ISHE, and it also means 
a possible application in the future. If we can take the 
direction of the spin Hall current as a new sign, i.e, +x di- 
rection means "0" and —x direction means "1" , and these 
two states maybe the basic for a new logic electronic- 
device 1 ^. 

The resonant ISHC predicted by usA 2 is caused by the 
energy crossing near the Fermi energy due to the tran- 
sition between mostly spin — i holes and spin-| holes. 
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FIG. 3: (a) The ISHC as a function of hole densities for 
B = 0.82T and am/fi 2 =0.05nm. (b) Resonant ISHC within 
Luttinger model including Zeeman splitting, and parameter 



are taken from Refi— except = 2.9 x 10 1 
g s = -0.44. 



and 



However, for the 2D heavy-hole system described by a 
fc-cubic Rashba model, we can deduce from Eq. (3) that 



(6) 



is steady for any n and Rashba coupling constant in the 
absence of Zeeman splitting. Then there is NO resonant 
ISHC within Hamiltonian (2) in the presence of a dc field. 

Besides, the resonant ISHC was also predicted in 2D 
quantum wells in a strong perpendicular magnetic field, 
where the Zeeman splitting plays a crucial role. The k- 
linear Rashba coupling, generated by spin-orbit interac- 
tion in wells lacking bulk inversion symmetry, competes 
with Zeeman splitting to introduce additional degenera- 
cies between different Landau levels (i.e., two states E nfl , 
E n ±i,±fi) at certain magnetic fields. This degeneracy, 
if occuring at the Fermi level, gives rise to a resonant 
ISHC 14 >. The Hamiltonian for a 2DEG in a perpendicu- 
lar magnetic field is 



He = 7Tl(P+—) 

lm c 



A. eA. 
-z ■ {p+ — ) 
h c 



x a 



(7) 



where m' , g' s are the electron's effective mass and Lande- 
g factor, respectively. In addition, /ib is the Bohr mag- 
neton, A is the Rashba coupling. The eigen energy of H e 
is given by 



= hio. 



(« + fv / (l-.9') 2 +8m 7 2 ) 



(8) 



where tu e — eB/m'c with fj, = ±1, for n > 1; and fi = 1 
for n = 0. Numerically, the degeneracies between differ- 
ent Landau levels {E n ^, £ n ±i,±/i) request that 



< g' < 1, 



(9) 



for a reasonable value of Rashba coupling X/h. To a turn, 
a set of parameters appropriate for Ino.53Gao.47As/In 
o.52Alo.48Asi2i are in this range for g' — 0.1. 

Let's discuss the effect of Zeeman splitting on the ISHC 
of a 2DHG. The 2D heavy-hole described by fc-cubic 
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Rashba model in a perpendicular magnetic field including 
Zeeman splitting can be written as 



H z = hio[(a'c 



1, 



■ / 3 4- t d 

-«7(a a —a' 



-g s ^bBa z (10) 



and the energy levels in units of fkv can be expressed as 



E ns = (n - 1) + - v /(3 - 3. 9 ) 2 + 4 7 2 n(n-l)(n-2), 

(11) 

where g = g s m/2m e . The key difference between fc-linear 
and fc-cubic Rashba model is the constant along with g, 
and "3" in (3 — 3<7) 2 is due to fc-cubic Rashba term, which 
is similar as "1" in (1 — g') 2 arising from fc-linear Rashba 
term. From Eq.(12), we can deduce that if the energy 
crossing between E ns and E n ±\^ s occurs, which may 
give rise to a resonant ISHC, requests that | < g < 1 
for a reasonable value of am/h 2 . In calculation, we take 
the effective mass m = 0.27r7i e (Rcf. 6 - 19 ), and then the 
corresponding effective gyromagnetic factor g s is in the 
range of 4.93 ~ 7.4. A typical value of g s has been given 
by Winkler—, and g s — —0.44. The requested g s is far 
away from the typical value of g s = —0.44, which means 
that the Zeeman splitting is too small to introduce "effec- 
tive" degeneracy between different Landau levels. So the 
resonance effect stems from energy crossing of different 
Landau levels near the Fermi level due to the competi- 
tion of Zeeman splitting and fc-cubic Rashba SOC is out 
of reach for real materials or in experiments. 

In addition, the Zeeman splitting plays a crucial role on 
the ISHC within fc-linear Rashba model. When Zeeman 
energy is neglected, the zero ISHC persists in a magnetic 
fieldSi. Within the perturbation method in this paper, 
we can also achieve this conclusion. However, the Zee- 
man splitting has much less effect on the ISHC of a 2DHG 
in magnetic field. To demonstrate this, a resonant ISHC 
within Luttinger model including linear Zeeman split- 
ting term has been shown in Fig. 3 (b), and parameter 
used are taken from Ref Ai£ except = 2.9 x 10 16 m~ 2 , 
and g s = —0.44. The resonance appears at B = 29.98T. 
Therefore, the resonant ISHC of a 2DHG in a perpendicu- 
lar magnetic field contributed from the interplay between 
mostly spin — i holes and spin-| holes is confirmed. 

To compare with the previous resuh— in the absence 
of magnetic field, we calculate the interband contribution 
to the ISHC as a function of the length scale am/h 2 of the 
Rashba coupling at a hole density of % = 3x 10 14 m~ 2 
for a fixed magnetic field. Our numerical result has been 
shown in Fig. 4. The circle red dots indicates v = 100 and 



the blue triangle indicate v = 1000. The ISHC in the ab- 
sence of magnetic field as a function of am/h 2 at the 
same hole density has been plotted in this fig too, which 
is denoted by the solid black line. To take a overview ex- 
pediently, it has also been shown in the fig inset solely and 
data are obtained within the method in RefA As seen in 
Fig. 4, the interband contribution to the ISHC increases 
with increasing Rashba coupling, and it approaches the 

value in absence of magnetic field in weak magnetic limit. 
15 




am/(h/2K) / (nm) 

FIG. 4: (Color online) The interband contribution to the 
ISHC as a function of am/h 2 at nn — 3 x 10 14 m -2 for a 
fixed magnetic field. The circle red dots indicates v — 100 
and the blue triangle indicates v — 1000. Solid line and inset: 
The ISHC as a function of am/h 2 at the same hole density in 
the absence of magnetic field, and data are achieved within 
the method in Ref^. 



In conclusion, two likely ways to establish ISHE in ex- 
periments shall be suggested. (1) A resonant ISHC shall 
appear in a 2DHG for a high hole density and in a strong 
perpendicular magnetic field, and this phenomenon may 
be used to distinguish the ISHE from the extrinsic one. 
(2) The ISHC shows rich sign changes, and a notable 
sign changes can be achieved for a lower hole density and 
moderate magnetic field. The sign of ESHE depends on 
the sign of scattering potentia l 1 ^ and such difference 
may be used to distinguish ESHE from intrinsic contri- 
butions and determine the intrinsic origin of the effect. 
The sign issue should be regarded as a very important 
aspect of SHE, and can be used in future experiments or 
application. 

We thank Shou-Cheng Zhang, Shun-Qing Shen and 
Wuming Liu for many helpful discussions. 
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